Preprint typeset in JHEP style - HYPER VERSION 



Anomalies and the chiral magnetic effect in the 
Sakai-Sugimoto model 



Anton Rebhan, Andreas Schmitt, Stefan A. Strieker 

Institut fiir Theoretische Physik, Technische Universitdt Wien, IO4O Vienna, Austria 



rebhcinaOhep . 


itp . 


tuwien. ac . 


at 


aschmittOhep 


. itp 


. tuwien.ac 


.at 


strickerOhep 


. itp 


. tuwien. ac 


.at 



Abstract: In the chiral magnetic effect an imbalance in the number of left- and right- 
handed quarks gives rise to an electromagnetic current parallel to the magnetic field pro- 
duced in noncentral heavy-ion collisions. The chiral imbalance may be induced by topo- 
logically nontrivial gluon configurations via the QCD axial anomaly, while the resulting 
electromagnetic current itself is a consequence of the QED anomaly. In the Sakai-Sugimoto 
model, which in a certain limit is dual to large- A'^c QCD, we discuss the proper implemen- 
tation of the QED axial anomaly, the (ambiguous) definition of chiral currents, and the 
calculation of the chiral magnetic effect. We show that this model correctly contains the 
so-called consistent anomaly, but requires the introduction of a (holographic) finite coun- 
terterm to yield the correct covariant anomaly. Introducing net chirality through an axial 
chemical potential, we find a nonvanishing vector current only before including this coun- 
terterm. This seems to imply the absence of the chiral magnetic effect in this model. On 
the other hand, for a conventional quark chemical potential and large magnetic field, which 
is of interest in the physics of compact stars, we obtain a nontrivial result for the axial 
current that is in agreement with previous calculations and known exact results for QCD. 
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1. Introduction 

Topologically charged gauge field configurations in QCD generate chirality due to the non- 
abelian axial anomaly. In the presence of a magnetic field, this chirality, i.e., an imbalance 
in the number of right- and left-handed quarks, has been predicted to generate an electro- 
magnetic current parallel to the applied magnetic field. This is a consequence of the QED 
axial anomaly and has been termed chiral magnetic effect ||l|, |2|, |3|. As a result, electric 
charge separation may occur in noncentral heavy-ion collisions, where magnetic fields up 
to 10^^ G can be generated temporarily, and corresponding experimental evidence has in 
fact been reported in refs. [Q, ^ (see however ref. |^). 

In a simplified picture, one may study the induced current for a static magnetic field. 
The generalization to time-dependent magnetic fields, as produced in heavy-ion collisions, 
in principle amounts to computing a frequency-dependent conductivity 0. However, the 
observed charge separation is proportional to the zero-frequency limit In this paper, the 
currents we compute always correspond to the zero- frequency limits of the conductivities. 
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Another simplification of the highly nontrivial scenario of a heavy-ion collision is to 
mimic the (event-by-event) topologically induced chirality by a nonzero axial chemical 
potential ^5, the difference of right- and left-handed chemical potentials. The resulting 
current is a vector current proportional to fi^. In a more general setup, although negligible 
in the heavy-ion context, one may also include a quark chemical potential ^u, which is the 
same for right- and left-handed fermions. Again via a nonzero magnetic field, an axial 
current is generated in this case ^. This effect may be of relevance for the physics 
of compact stars [^], where strongly interacting matter can reach densities of several 
times nuclear ground state density, and (surface) magnetic fields up to 10^^ G have been 
measured, indicating the possibility of even higher magnetic fields in the interior. Also the 
direct high-density analogue of the chiral magnetic effect has been studied in the context 
of neutron star physics |11|. 

In the present paper, we apply a strong-coupling approach, based on the AdS/CFT 



correspondence [l^, |13|, to compute both kinds of currents. We use a general setup to 
account for nonzero temperatures, relevant in the context of heavy-ion collisions, as well 
as for nonzero quark chemical potentials, relevant in the astrophysical context. Besides 
the chirally symmetric phase we also consider the chirally broken phase which is important 
in both contexts: heavy-ion collisions are expected to probe the region of the QCD chiral 
phase transition; in quark matter at densities present in compact stars, chiral symmetry 
may also be spontaneously broken, for example in the color-flavor locked phase [|l5| ]. 

We use the Sakai-Sugimoto model fl^, which is sometimes called "holographic 
QCD" since in the limit of small 't Hooft coupling it provides a string theory dual to 
large- A'^c QCD. However, we work in the opposite, strongly coupled, limit where the sim- 
ple gravity approximation can be employed but where the model is no longer dual to 
QCD. The model still yields interesting qualitative predictions especially in view of the 
strong-coupling nature of both contexts mentioned above, i.e., in QCD at large (but not 
asymptotically large) temperature and small quark chemical potential and QCD at large 
(but not asymptotically large) quark chemical potential and small temperature. 

The Sakai-Sugimoto model is particularly suited for our purpose since it has a well- 
defined concept for chirality and the chiral phase transition. It is straightforward to in- 
troduce right- and left-handed chemical potentials independently. Several previous works 
have considered currents in a magnetic field at nonzero chemical potentials in this model 
1 18, 19, 20, The purpose of the present paper is two- fold. The physical motivation is 
to extend these calculations to the currents relevant for the chiral magnetic effect, and to 
compare our strong-coupling results to the weak-coupling results as well as the lattice 
results [22 1 in the existing literature. There is also a more theoretical purpose of our work, 
addressing certain fundamental properties of the Sakai-Sugimoto model. We discuss in 
detail how to implement the covariant QED anomaly into the model in order to obtain 
physically acceptable predictions. Moreover, we elaborate on an ambiguity in the defini- 
tion of the chiral currents in the presence of a Chern-Simons term that has been observed 
previously H, H, ^ (see also @). 

Our paper is organized as follows. We start with a brief introduction into the model 
and a general discussion of the currents, in particular the appearance of consistent and 
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covariant anomalies, in sec. ||. In sec. ^ we discuss the solution of the equations of motion 
in the presence of background magnetic and electric fields. We present analytical solutions 



for the chirally broken phase, sec. 3.1, and the chirally symmetric phase, sec. |3.2| . We then 
discuss the ambiguity of the currents, defined on the one hand via the general definition 
from sec. ^, and on the other hand from the thermodynamic potentials obtained in sec. |3|. 
In sec. ^ we present our results for the axial and vector currents and give our conclusions 
in sec. ^. 



2. Anomalies in the Sakai-Sugimoto model 
2.1 Brief introduction into the model 

The Sakai-Sugimoto model is based on ten-dimensional type-IIA string theory, with a 
background geometry given by Nc D4-branes. They span four-dimensional space-time (r, x) 
and a fifth extra dimension X4 compactified on a circle whose circumference is parametrized 
by the Kaluza-Klein mass Mkk; 2:4 = X4 -|- 27r/MKK- Through this compactified dimension 
and antisymmetric boundary conditions for fermions supersymmetry is completely broken. 
Left- and right-handed chiral fermions are introduced by adding Nf D8- and Nf D8-branes 
which extend in all dimensions except X4. In this compact direction, they are separated by 
a distance L € [0, tt/Mkk]- For more details about the setup of the model and the explicit 
form of the background metric, including the holographic direction z and a four-sphere 5^, 
see the original papers by Sakai and Sugimoto, refs. [p!6| , 0]. We employ the probe brane 
approximation, i.e., the background geometry is assumed to be unaltered by the flavor 
branes. This is a good approximation for Nf <^ Nc- 

There are two possible background geometries. One, interpreted as the confined phase, 
has a cigar-shaped {xi,z) subspace, ending at the tip z = 0, and a cylinder-shaped (r, z) 
subspace, where r is Euclidean time on a circle with circumference given by the inverse 
temperature, t = t + 1/T. In the other geometry, interpreted as the deconfined phase, 
X4 and r exchange their roles, such that the {x/^,z) subspace is cylinder-shaped while the 
(r, z) subspace is cigar-shaped, corresponding to a geometry with black-hole horizon after 
analytical continuation ir ^ t. 

Chiral symmetry breaking is realized in the model as follows. A U (Nf) gauge symmetry 
on the fiavor branes corresponds to a global U{Nf) at the boundary. Therefore, the bulk 
gauge symmetries on the D8- and D8-branes can be interpreted as left- and right-handed 
flavor symmetry groups in the dual fleld theory. The Chern-Simons term accounts for the 
axial anomaly of QCD, such that one is left with the chiral group SU{Nf)L x SU{Nf)ji 
and the vector part C/(l)y. There is no explicit breaking of this group since the model 
only contains massless quarks. Spontaneous chiral symmetry breaking is realized when 
the D8- and D8-branes connect in the bulk. They always connect in the confined phase, 
where the {x4,z) subspace is singly connected. Whether they connect in the deconfined 
phase depends on the separation L of the D8- and D8-branes in the extra dimension X4. 
Here we shall always consider maximally separated branes, L = tt/Mkk- With this choice 
the fiavor branes necessarily extend to the black-hole horizon and thus never connect in 
the deconfined phase. Consequently, the deconfinement and chiral phase transitions are 
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identical and happen at a critical temperature Tc = MKK/(2vr), i.e., when the radii of 
the r and X4 circles are equal. We shall always use the terminology of chirally symmetric 
and chirally broken phases. This is equivalent to speaking about deconfined and confined 
phases, but more appropriate in our context because we are interested in the interplay of a 
magnetic field with chirality. In the probe brane approximation, the chiral/deconfinement 
phase transition is given solely by the background geometry and is not affected by the gauge 
fields on the flavor branes. In particular, it does not depend on the chemical potential, 
which, at least for the deconfinement phase transition, is in accordance with expectations 
for large- A?^, QCD ||]. 

2.2 Action, equations of motion, and currents 

In this section we discuss the general equations of the model in the broken phase where the 
D8- and D8-branes are connected. The equations for the symmetric phase are very similar 
and shall be given later where necessary. The D-brane action consists of a Dirac-Born-Infeld 
(DBI) and a Chern-Simons (CS) part. We approximate the DBI action by the Yang-Mills 
(YM) action which is a good approximation for small magnetic fields. The use of the YM 
action greatly simplifies the treatment since the equations of motion then have solutions 
which can be given in an almost entirely analytical way. Throughout the paper we shall 
work with one quark flavor, Nj = 1. The currents we compute are expected to be simple 
sums over quark flavors, each flavor contributing in the same way, distinguished only by 
its electric charge. This is rather obvious in the chirally symmetric phase. In the chirally 
broken phase, the flavor contributions may be more complicated in the case of charged 
pion condensation. However, since we work at vanishing isospin chemical potential, there 
is only neutral pion condensation and the different flavor contributions decouple. 
For one quark flavor and the gauge Az = the (Euclidean) action 



S = 5ym + "^cs (2.1) 



is given by ll 



•S'ym = /^-^KK / I dz 



00 



00 

00 



(2.2a) 



5cs = d^x j dz Af,Fz,Fp^e^'''P% (2.2b) 

with Greek indices running over /i, i/, . . . = 0, 1, 2, 3. Our convention for the epsilon tensor 
is £0123 = +1- In eq. ( |2.2| ) we have defined the metric functions 

k{z) = l + z\ h{z) = {1 + z^)-^/^ , (2.3) 

and the dimensionless constant 

(2.4) 



2167r3 ' 

where A is the 't Hooft coupling. The integration over the four-sphere has already been 
done, and we are left with the integral over space-time (r, x) and the holographic coordinate 
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z which extends from the left-handed boundary (z = -|-cxd) over the tip of the cigar-shaped 
(x'^jz) subspace (z = 0) to the right-handed boundary (z = — oo). The coordinate z is 
dimensionless and is obtained from the dimensionful coordinate z of ref. upon defining 
z' = z/ukk {T < Tc) and then dropping the prime. Here, ukk = 4i?^M^j^/9, with R 
being the curvature radius of the background metric. Since we work at finite temperature, 
we need to work in Euchdean space. However, we use Minkowski notation which is more 
convenient for the following discussion of the anomaly. More precisely, we start from the 
Euclidean action with imaginary time r and replace ~^ ^^O) after which we may write 
the result using a Minkowski metric with signature (— , -|-, -|-, -|-). The space-time integral 
is denoted by d'^x for simplicity but actually is an integral dr d^x over imaginary time r 
and three-dimensional space. For the general form of the action, without any gauge choice 
and for more flavors, see for instance refs. |1^, 2C, The equations of motion for Nf = 1 
are 



KMl^d-,[k{z)F''^'] + Kh{z)duF''f' 
KM^j^d^[k{z)F'''] 



Nr 



_Nc 
647r2 



(2.5a) 
(2.5b) 



where the second equation is obtained from varying Az in the action prior to setting = 

Next we introduce the chiral currents. The usual way is to define them through the 
variation of the on-shell action with respect to the boundary values of the gauge fields (see 
however ref. ||2^ for a discussion of possible alternatives). We thus replace A^{x,z) — > 



A^(x^ z) + (5A^(x, z) in the action and keep the terms linear in 5A^^ 



[x,z 



to obtain 



6Sym = 2/tM^K i / d^xkiz)F^f'6A^ 



+ I d'x I dz^F-^^dA^ 

-co ^"^KK 



d X I dz 

'OO 

Nr ( 1 



h{z) 



d,[k{z)F^'^] + ^d,F'^'^ 



6A,, 



SScs = ^\-^ I d'^xA.FpJA^ 



KK 

2 r roo 

-jd^xj dzA^FzuSA^ 



+ I d^x I dzFz,Fp^6Af,}e''''P'' 



(2.6a) 



(2.6b) 



In the total variation 5S = 5Sym + '^'S'cs the bulk terms vanish upon using the equation 
of motion for Afj_ ( |2.5a ) and we are left with boundary terms only. According to the 
holographic correspondence, we keep only the boundary terms at \z\ = oo and drop any 
terms from space-time infinities. This may seem natural but possibly is problematic in our 
case as we shall discuss later after we have implemented our specific ansatz. The boundary 
terms at the holographic boundary z = ±00 lead to the left- and right-handed currents 



q-p 



5S 



6Afj_{x, z = ±00) 



(^2KMl^kiz)F^'^ - -^eM-PM^Fp.) 



(2.7) 



2 = ±00 



where the first (second) term is the YM (CS) contribution. This result of the currents is in 
agreement with refs. [pO| , [2^ , 26 1, see also p7| , [2^ . The overall minus sign in the definition 
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originates from our use of the Euclidean action which is minus the Minkowski action, and 
the functional derivative is taken with respect to the space-time coordinates x (and not 
also with respect to the holographic coordinate z plus a subsequent limit z — > ±00). The 
currents (|2.7D can also be obtained from 



'L/R QdA 



(2.8) 

2 = lt0O 



in accordance with the usual rules of the gauge/gravity correspondence. 

As already pointed out in ref. it is only the YM part of the current, i.e., the first 
term in eq. (^), which appears in the asymptotic expansion of the gauge fields. From the 
definition ( p.T] ) and with k{z) = 1 + we find 

A,(.,.) = A„(.,. = ±oo)±^- + o(^). (2.9) 

One can also confirm this relation from our explicit results in the subsequent sections. 
2.3 Consistent and covariant anomalies 



The divergence of the currents (|2.7| ) can be easily computed with the help of the equation 
of motion for ( 2.5b| ). One obtains 



d,J^/R = MJym + JcsTl/u = (1 - I) F^J^'n/R ' (2-10) 

with the left- and right-handed field strengths F^,) (x) = F^y{x, z = ±00), and the left- and 
right-handed dual field strength tensors = \ FpJ^e^'^^" . (For notational convenience 
we use the labels L, R and related labels such as V, A sometimes as superscript, sometimes 
as subscript.) With the vector and axial currents 

j'^ ^ + , ji; ^j^- , (2.11) 

and the vector and axial field strengths introduced as = -|- F^, F^^, = — F^, 
eq. ( p.lOP yields the vector and axial anomalies 

d,J' = -^F^.Fa , (2-12a) 
« = ^ (<^r + ^A^r) • (2.12b) 

The coefficients on the right-hand side (which as we saw receive contributions from both the 
YM and CS parts of the currents) are in accordance with the standard field theoretic results 
for Nc chiral fermionic degrees of freedom coupled to left and right chiral gauge fields | |29| . 
The above form of the anomaly, which is symmetric in vector and axial- vector gauge fields, 
is called consistent anomaly. If left- and right-handed Weyl spinors are treated separately, 
this form of the anomaly arises unambiguously. This is explained for instance in ref. |^0| , 
where left- and right-handed fields are separated by an extra dimension. This is not unlike 
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our present model and it is thus not surprising that the consistent anomaly arises naturally 
from the above definition of the currents. In QED, however, we must require that the 
vector current be strictly conserved, even in the presence of axial field strengths. As was 
first discussed by Bardeen [^], this can be achieved by the introduction of a counterterm 
that mixes left- and right-handed gauge fields. Having even parity, Bardeen's counterterm 
is uniquely given by |30| 



AS = c I d^A^lA^Fp. + A^^A^F^„)e^''f"^ , (2.13) 

where c is a constant determined by requiring a strictly conserved vector current. Because 
this expression can be naturally written as a (metric-independent) integral over a hyper- 
surface at l^l = A ^ oo with left- and right-handed fields concentrated at the respective 
brane locations, A.S can actually be interpreted as a (finite) counterterm in holographic 
renormalization. In particular, it does not change the equations of motion. 

To obtain the contribution of Bardeen's counterterm to the chiral currents we replace 
a]i'^ ^ a]1^^ + M^/^ to obtain 

5^S = ±cj d^x [a^/'^F^J'^ - A^I^'Fl^J'^ + 2A«/^F^V'') <5^^/^6^"'"^ 



2c / d^a^^^^^V^^^VR 



gM^'P'^ . (2.14) 



Again dropping the space-time surface terms, the contribution to the currents is therefore 
^^L/R = [a^^'-F^J'^ - A^/'^F^J'^ + 2A^/^F^V^) e^^/"^ , (2.15) 
and the contribution to the divergence of the currents becomes 

= Tc {F^TFi^/L + F,'J^Fi^/L) ■ (2-16) 
Denoting renormalized left- and right-handed currents as 

•^L/R = "^L/R + ^^LjR ' (2-17) 

and similarly the renormalized axial and vector currents as ^7^, we find that the choice 



487r2 

leads to the covariant anomaly 



(2.18) 



df,J'' = 0, (2.19a) 



Note that the prefactor in front of the first term in the axial anomaly now has changed 
to iVc/(87r2), from iVc/(247r2) in eq. ( p.l2b| ), which is the well-known result for the Adler- 
Bell-Jackiw anomaly for QED fsTI, |32] and which is essential for getting the correct pion 
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decay rate vr*^ 27. The necessity of adding the counterterm ( 2.13 ) to the Sakai-Sugimoto 
model is in fact completely analogous to the very same and well-known procedure in chiral 
models where a Wess-Zumino-Witten term accounts for the anomaly ||3^ . 

In the literature sometimes the coefficient of the subleading term in the asymptotic 
behavior of \z\ 00) and thus the YM part of the current (see eq. ( |2.9D ) is identified 

with the full current Q, see also [34, 35, Using this identification, it has also been 



assumed that the equation of motion for Az ( p. 5b ) represents the anomaly equation |21 



Indeed, from eq. ( 2.5bD one obtains the apparent anomaly 

'ym,l/r ''iqtt'^ i^" l/r 

which leads to 



B 7^ - mJ^F^/^F'^" (2 20) 



d^^YM,, = ^ (<.^r + ^A^r) , (2.21b) 

and this does contain the same coefficient in front of F^jjFy' as the full covariant anomaly 
(2.19). However, it differs from the latter in the presence of axial gauge fields. In particular. 



the vector current is then not strictly conserved. The renormalized current Jl/r satisfies 
eq. ( 2.20 ) only for = Fj^^, } Even when this issue may be ignored, because all axial 



vector field strengths are set to zero, it appears to be questionable to keep only part of the 
full current (^). 

In the remainder of the paper we shall consider the full currents for which Bardeen's 
counterterm is needed, and study the implications, which indeed differ from keeping only 
the YM part of the currents. (The effect of truncating to the YM part can be easily read 
off from the expressions that we shall give.) 

3. Background electromagnetic fields and chemical potentials 

The discussion in the previous section was general in the sense that we have not specified 
any gauge fields except for the gauge choice = 0. In this section we specify our ansatz ac- 
cording to the physical situation we are interested in. This includes a background magnetic 
field B as well as separate left- and right-handed chemical potentials /iL,_R = Aq{z = ±00), 
or, equivalently, the ordinary quark chemical potential ^ = {hr + iil)/^ and an axial chem- 
ical potential ^5 = (/i/j — fJ-i)/'^- With these ingredients we can obtain results relevant for 
the heavy-ion context (nonzero ^5, negligibly small /u) and for the astrophysical context 
(vanishing //5, large fi). In order to be able to check the axial anomaly explicitly, we also 
add an electric field E and an "axial electric field" e parallel to the magnetic field. The 
electric field E is needed because the axial anomaly is proportional to E • B. The (unphys- 
ical) field e shall be used to check the absence of a vector anomaly, i.e., the conservation 



The more general validity of eq. (2.20) has been assumed incorr ectly in eq. (2.1) of ref . and eq. (36) 
of ref. i. 
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of the vector current, which must be true even in the presence of e, and would be trivial 
without e. In our final results for the currents, the electric fields are however set to zero. 

For previous discussions of background electric and magnetic fields in the Sakai- 
Sugimoto model see for instance refs. 2C, We shall only consider spatially 



homogeneous systems. This is the simplest case, which might however require generaliza- 
tion when the true ground state is more complicated, for instance when Skyrme crystals 



are formed |4C] 



3.1 Chirally broken phase 

In our ansatz the nonzero fields are ^o(^) z), Ai{x2), ^sit, z), where the dependence on t will 
only be present for nonvanishing electric fields E and e at the holographic boundary. The 
temporal component is needed to account for nonzero (left- and right-handed) chemical 
potentials which correspond to the values of Aq at the boundary. The electromagnetic fields 
are encoded in the boundary values of the spatial components. Since the gauge symmetry 
in the bulk corresponds to a global symmetry for the dual field theory, the fields at the 
holographic boundary are not dynamical and merely serve as background fields. This is 
however sufficient for our purpose. The magnetic field B is assumed to point into the 
3-direction, B = (0,0, -B). Consequently, we can choose 

Ai{x2) = -X2B (3.1) 

at the holographic boundary. The equations of motion show that Ai can be chosen to be 
constant in z throughout the bulk. (This is different in the presence of an isospin chemical 



potential [^.) Consequently, Fi2{x, z) = B. For notational convenience we have absorbed 
the electric quark charge qj into B, i.e., actually B qjB with qj = 2/3 e for f = u, and 
qf = — l/3e for f = d. With nonzero Aq and Ai, accounting for the chemical potential 
and the magnetic field, a nonzero ^3 is induced, even without electric field. In the broken 
phase, ^3 develops a nonzero boundary value, corresponding to the gradient of the neutral 
pion |2^. Just as for a usual superfluid, where the gradient of the phase of the 

order parameter is proportional to the superfluid velocity, this gradient of the pion field 
can be viewed as an axial super current pO|]. 

Next we introduce the electric field E = (0, 0, E) parallel to B and, as explained above, 
an "axial electric field" e = (0, 0, e). We thus have to add —t{E^ e) to the boundary value 
of A3, such that, together with the axial supercurrent jt, we have 

A3{t,z = ±^) = -t{E^e)Tjf (3.2) 

Due to the axial electric field we allow the supercurrent to become time-dependent. Strictly 
speaking the electric fields prevent us from using a thermodynamic description since it 
introduces a time-dependence and thus non-equilibrium physics. Therefore, our electric 
field should be considered infinitesimal. This is sufficient for our purpose since we can 
check the anomaly relations with an arbitrarily small electric field. Moreover, as mentioned 
above, the physical situations we are interested in do not require finite electric fields anyway. 



-9- 



With the above ansatz, the YM and CS contributions to the action ( |2.2| ) become 

(3.3a) 



/f'OO 
<fx dzk[z)[-{d,AQf + [d,A^) 
J —CO 



Scs 



12^2 



dz 



- Ai [{d2A^){d,A^) - {d2A^){d,A^)] } . (3.3b) 



We have written all terms which are needed to derive the equations of motion, including 
the ones that vanish on-shell. More specifically, the second line in the CS action ( |3.3b ) 
vanishes on-shell because neither nor ^3 depends on X2, but yields a finite contribution 
to the equations of motion. The equations of motion are 



d,{kd,A^) = 2pd,Ao, 

dtikdM = 2(3dtA3, 



with the dimensionless magnetic field 



aB 



(3.4a) 
(3.4b) 
(3.4c) 



(3.5) 



and a = 277r/(2A). We defer the details of solving the equations of motion to appendix 
)r the g 

Ao{t,z) = fit- IJ'5,t- 



A.l. The results for the gauge fields and field strengths are 

sinh(2/3 arctan z) 



A3it,z) 

and 

kd.Ao = -2(5 
kd.As = -2(5 
Here we have denoted 



-{Jt-et) 
-tE - fi^^t 
-ijt-et) 



sinh (5tt 
cosh (2/? arctan z) 
sinh /Jtt 
cosh (2/3 arctan z) 
sinh Ptt 
sinh(2/3 arctan z) 



coth Ptt 
- coth p-K 



sinh Pn 
cosh (2/3 arctan z 



sinh Pn 
sinh (2/3 arctan z) 
sinh p-TT 



+ {Jt - et) 
+ {jt - et) 



sinh (2/3 arctan z) 

sinh Pn 
cosh(2/3 arctan z) 

sinh p-TT 



fit = 1^ + et coth Pit , fi^^t = fi5 + Et tanh Pn . 



(3.6a) 
(3.6b) 

(3.7a) 
(3.7b) 

(3.8) 



i.e., both boundary values of Aq become time-dependent through the electric fields. As can 
be seen from the detailed derivation in appendix |A.l , this time-dependence is unavoidable. 
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Figure 1: Gauge fields in the chirally broken phase as functions of the holographic coordinate 
z e [— cxDjCxd] for a finite quark chemical potential and vanishing axial chemical potential (left) 
and vice versa (right). Dashed lines: gauge fields with vanishing magnetic field; solid lines: gauge 
fields with a nonzero magnetic field /3 = 0.6. In both plots we have set the electric fields to zero, 
= e = 0. The boundary values at 2; = ±00 correspond to left- and right-handed quantities. The 
magnetic field induces an axial supercurrent (boundary value of ^3) in the case of a nonvanishing 
quark chemical potential. If both [i and /^s are nonvanishing, the gauge fields are neither symmetric 



nor antisymmetric in z. The analytic expressions for these curves are given in eqs. (3.6) 



In fig. H we plot the gauge fields for = e = at the minimum, i.e., after jf has been 
determined to minimize the free energy, see below. 

The thermodynamic potential S7 = ^5'on-sheii is obtained from eqs. ( A.13| ), treating t 
as an external parameter. 



8kM? 



KK 



{ [ijt - etf - piP) + p [fitijt - et) + t^^s,*] } , (3.9) 



3 

where we have abbreviated 

7r/?2 f — + for p^O 

2^^^^^^ i \P\ for|/3|^oo 
Minimization of O with respect to jt yields the axial supercurrent 



Pp 



P coth Pit 

MP) . 



(3.11) 



We see that the supercurrent depends neither on p^ nor on E. Therefore, at t = it is 
simply the one-flavor limit of the result obtained in ref. ||2^ (where the D8 and D8 branes 
were identified with R and L, not with L and R, respectively, hence the different sign of 
the supercurrent). 

3.2 Chirally symmetric phase 



As explained in sec. 2.1, in the chirally symmetric phase the D8 and D8-branes are not 



connected. On both branes the holographic coordinate z now runs from z = 0, the black 
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hole horizon, to the holographic boundary at z = oo, and both branes yield separate 
contributions to the action, 



S — {Syu + ^Yu) + (•S'cs ~ ^cs) ■ 



(3.12) 



The CS action assumes different overall signs on the D8- and D8-branes since its parity is 
odd. The YM and CS contributions are 



= ^^KK^" / d^'x 



dz 



127r2 



ko{z){d,A'^f + h{z){d,Alf 
dz{(a24)[4(5,4)- 4(5,4)] 

4[(524)(5.4) - (924)(5.4)]} , 



with h = L,R. Here, we have defined the dimensionless temperature 

2ttT 



Mkk 



(3.13a) 



(3.13b) 



(3.14) 



In contrast to the broken phase there are different metric functions for temporal and spatial 
components of the gauge fields, 

(1 + ^2)3/2 



ko{z) 



ksiz) 



z(l + z 



2a1/2 



(3.15) 



Note the slight difference in notation of the gauge fields: while in the broken phase 
A^ (x) = Af^{x,z = iLcxd) always implies evaluation at the holographic boundary, here 

r / D 

we label the bulk gauge fields A^ [x, z) by L and R to indicate whether they live on the 
D8- or on the D8-brane. Since we always discuss broken and symmetric phases separately, 
this should not cause any confusion. 

The equations of motion on the separate branes become 

2/3, 



d.ikod.Al"') 
5.(^39,4/^) 
dtikod.A^J'') 



+—rl A^^^ 

2P L/R 
± ^3 "^^^O ' 



2/?, 



±^dtA^ 



L/R 



(3.16a) 
(3.16b) 
(3.16c) 



Details of solving the equations of motion are presented in appendix |A.2| . The final solution 
for the gauge fields is 

Po 



^o^^it^z) = {jit T /^5,t) 



p{z) 



■Qiz) 



4/^(t,z) 



kodzP - — (1 + kodzq) 
Qo 



(3.17a) 
(3.17b) 



which is plotted in fig. ^ for E = e = 0. Below we shall also need the field strengths on the 
branes. 



kodzA^^^ = {fit =F ^5,i) ( kodzp - —kod^q ) , 

V qo J 



p{z) 



Po 

— ( 

qo 



-q{z) 



(3.18a) 
(3.18b) 
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Figure 2: Left- and right-handed gauge fields (left and right panel, respectively) in the chirally 
symmetric phase as functions of the holographic coordinate z £ [Q,oo] for /i — 0.9, 115 — 0.1. We 
have set the electric fields to zero, E = e = 0. The temporal components Aq^^ approach the 
chemical potentials ^^115 at the boundary z — 00, while the spatial components A^^^ vanish at 
z = 00. A finite magnetic field (solid lines, here (3/6^ = 0.6) distorts the gauge fields compared to 
the case of vanishing magnetic fields (dashed lines). In particular, the spatial component develops 
a nonzero value at 2: = 0. The different sign of this value for left- and right-handed fields, i.e., on 
the D8- and D8-branes, ensures the correct parity behavior of the fields. The analytical expressions 
for these curves are given in eqs. ( |3.17| ). 



The functions p{z), q{z) are hypergeometric functions which we defined in eqs. (A. 19) and 
which depend on the ratio /3/0^. Their values at z = are denoted by pq, qo, see eqs. 
( A.22| ), and the ratio po/qo behaves for small and large magnetic fields as 

' 1 + (2/3/^3)2(1^^4 for/?/03 

2\p\/e^ for \p\/e^ - 



Po 

— — < 



(3.19) 



00 



The boundary values of the temporal components are A^^^it, z = 00) = /J-tT /^5,t with 

fit = f^ + 2te—— , fj,5^t = fJ-b + ^tE—— . (3.20) 
Po d-^ Po 

It is instructive to compare this behavior of the axial chemical potential with the expected 
behavior for free fermions in a magnetic field. To this end, consider the lowest Landau 
level in which the spin of all (say, positively charged) fermions is aligned parallel to the 
magnetic field. As a consequence, all right- (left-) handed massless fermions move parallel 
(antiparallel) to the magnetic field. An electric field parallel to the magnetic field now 
shifts all momenta in the positive 3-direction by an amount Et. Consequently, some of the 
left- handed fermions are converted into right-handed fermions and a shift Et is induced 
in the difference of right- and left-handed Fermi momenta, (p^ — p^)/2 = Et |2|, 
Interpreting /is^t as (p^— p^)/2 (strictly speaking there is no well-defined Fermi momentum 
in our model), eq. (|3.20 ) reproduces this shift for asymptotically large magnetic fields 



because in this case qo/po ^^/(2/3)- For small magnetic fields qo/po — >■ 1, and the shift 
becomes linear in the magnetic field. Since oc B/T^, we can in principle also obtain 
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A*5,t = 1^-5 + tE for sufficiently small temperatures and fixed magnetic field. However, 
we cannot reduce the temperature arbitrarily in the above expression since below the 
critical temperature Tc we are in the chirally broken phase. In this case the analogous, 
temperature-independent relation in eq. ( ^ ) holds. 

The free energy, obtained from the YM and CS contributions ( |A.33 ), is 



(3.21) 



where we introduced the function 



r,{p/e^)^Io-{2P/d'fh + 2^ 

Qo 



3 + {2/3/e^y{ln4-l) for/3/^3^0 



for \P\/e^ oo 



(3.22) 



with integrals Iq and /s defined in eqs. ( |A.34 ) 



3.3 Ambiguity of currents 

In the following discussion we restrict ourselves to the symmetric phase, but one can easily 
check that all arguments hold for the broken phase as well. Let us first give the analogue 
of the definition of the currents (|2.7D for the symmetric phase. 



(3.23) 



where the notation fcj-^) (no summation over ^u) indicates the different metric functions for 



temporal and spatial components, see eq. ( p. 15 ). Equivalently, and in analogy to eq. (^), 
we can write the currents in the symmetric phase as 



dC 



(3.24) 



We shall now show that the currents defined via these equations are different from the 



ones obtained via taking the derivative of the free energy (3.21) with respect to the corre- 
sponding source. We do so for the vector density, i.e., the sum of left- and right-handed 
0-components of the currents. One can observe the same ambiguity for the other nonva- 
nishing components. The following arguments do not depend on the electric fields, so we 
temporarily set e = E = ^ for simplicity (and for a truly equilibrated situation). From the 
definition ( 3.23| ) and the gauge fields ( 3.17 ) and field strengths ( 3.1§| ) we obtain 



90 



(3.25) 



On the other hand, the free energy 17 of the system should yield the number density via 
the thermodynamic relation 



n 



9/i 



-firj. 



(3.26) 
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This result shows that n ^ which, given spatial homogeneity, is inconsistent. This 
inconsistency is absent for vanishing magnetic fields: using the behavior of the functions 
Po/lo and 77 from eqs. ( p.l9D and ( |3.22| ) one sees that for /? = the expressions for and 
n are identical. We can formulate this observation in a more general way. To this end we 
write the left- and right-handed on-shell Lagrangians, i.e., the integrands of the on-shell 
action ( 3.12| ), as ChiA^, d^A^, A^, dzA'i^), where all arguments of £/j depend on the chemical 
potentials with h = L,R and fJ-i/R = I^T /^5- Then, with = T/V f d^x dz Ch we 
have 



dfJ-h 



T 
V 



E 

i=0,3 



d^X 



dz 



dd,A^ 



T 



^ j=0,3 



d^X 



dCh dA 



2=0 



+ / d^x 



dz do 



dCh dAf 
dd2A'l dfih 



(3.27) 



where we have used partial integration and added and subtracted the derivative term in 
X2 in order to make use of the equations of motion. Now we use 



dAi^ 



dA^ 



dAl^ 



dd-,A\ 



0, 



(3.28) 



2 = 



2=0 ^i^h 

which follows from the explicit solutions (|3.171 ), whose behavior at 2; = 0, z = cxd is obtained 
with the help of eqs. (^]2l|), ([A^ ), and (^]2|). With these relat ions and the definition of 
the current from eq. ( 3.24[ ) we obtain 

dA^ 



d'^x I dzd2 ^^ ^ ; ■ (3.29) 

h=L,Ri=0,3 

This is the general form of the difference between the density defined as the 0-component of 
the current defined via eq. (3.24) and the density defined via the thermodynamic relation 
(3.26). For an explicit check of this relation one inserts the expressions 



82 



dC 



L/R 



KK 



L/R 



dC 



L/R 



dd2A^^^ 



and 



dA, 



L/R 



dl^L/R 

dA^^^ 



df^L/R 

into eq. ( 3.29| ). This yields 



p{z) 



2/3 



Po 



Pd,A 



Po 



kodzP (1 + ^0^2 

Qo 



dji 3 



wy J PO 
^) ''~7o 



(3.30a) 
(3.30b) 

(3.31a) 
(3.31b) 

(3.32) 
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With the definition ( ^.22 ) this confirms the difference between n and J^^ obtained from 
eqs. and (|]2|). 

From the general form ( p.29 ) we see that the additional term is a boundary term at 
the spatial boundary of the system. This suggests that the ambiguity in the currents is 
related to the terms we have dropped in sec. see eqs. (^.61). These terms correspond to 
currents at the spatial boundary and disappear in the presence of a homogeneous magnetic 
field only if the variation 5A^{x,z = ±00) can be chosen to vanish at this boundary. So 
this problem might be resolved by considering more complicated, spatially inhomogeneous 
gauge fields. In our homogeneous ansatz, it is however a priori not clear which definition 
of the currents corresponds to the correct physics. 



A possible solution to this ambiguity was suggested and applied in refs. [18, 21, 42 
In these references, the CS action has been modified according to 



Ih 

CS 



AT r fco f o p 

^ J d'x dz [4(5.4) -4(9.4) 



1 



(9,4) 4(^24) - 4(^24: 



(3.33) 



This modified action (marked by a prime) is obtained from the original CS action ( 3.13b| ) 
by adding a boundary term at the holographic and the spatial boundary, 



Q/h 



5h I 
CS "T '-'1 



boundary ' 



(3.34) 



with 



5, 



boundary 247^2 



dzA\ 



d^x4[4(a24)- 4(^24) 



z=oo 
z=0 



(3.35) 



Note that this boundary term cannot be considered as a holographic counterterm since 
it involves an integration over z. From eq. ( 3.33| ) we see that the addition of •S'boundary 
effectively amounts to a multiplication of the on-shell action by 3/2 because the second 
line in eq. ( |3.33| ) vanishes on-shell. The benefit of the modified action is that the integrand 



on the right-hand side of eq. (|3.2S| ) vanishes now, i.e., there is no ambiguity in the currents 
anymore. 

Modifications of a CS action by boundary terms are in fact sometimes necessary in 
order to ensure validity of the variational principle in the presence of nontrivial boundary 
values However, this is not what the above modification is achieving. Instead, it leads 
to gauge invariance under the residual gauge transformations Ai ^ Ai + 9iA(xi) which 
are compatible with the boundary conditions of our ansatz and which do not vanish at 
spatial infinity X2 = ±00 ||l^. In fact, by this modification one loses all anomalies for the 
(now uniquely defined) currents, as we show now. To this end, we switch on the electric 
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fields again. Then, the currents of the original action in the symmetric phase are 



T>Tr.: 



L/R 



J. 



2kM|k — 
0, 



L/R 



(3.36a) 
(3.36b) 
(3.36c) 
(3.36d) 



where we have used the definition (|3.23| ) and the gauge fields ( 3.17| ) and field strengths 



( 3.18 ). All terms containing a 1/3 originate from the CS contribution of the current, i.e., 
from the second term in eq. ( |3.23D . All other terms are YM contributions. In particular, 
the 2-component of the current is a pure CS term. This component is unphysical because it 
depends on our choice to introduce the magnetic field via the gauge field Ai. We could have 
introduced the same magnetic field via A2 or a combination of Ai and A2, in which case 
the 1- and 2-components of the currents would have been different. We shall see below that 
Bardeen's counterterm solves this problem by canceling the 2-component. Here, however, 
it gives a nonzero contribution to the anomaly. Namely, the divergence of the (unmodified) 
currents becomes 



Nr 



(3.37) 



where we have used the definition of i^^^t ( P-20|) and kM^-^P = NcB/Wtt"^. This is exactly 
the consistent anomaly (|2.1[l|) , because 



=F 



127r2 ^ ^ ' 



(3.38) 



The new currents iT^'^^ from the modified action are simply obtained by multiplying the 
CS contribution of the currents ( 3.36 ) by 3/2. Doing so in the explicit results ( |3.36| ), this 
yields 



/2 



(3.39) 



which can also be inferred in generality from (l2.1C| ). Consequently, the anomaly has dis- 
appeared. In other words, the new vector and the axial currents are both conserved. 
Nevertheless, one finds nonzero currents in the direction of the magnetic field. Multiply- 
ing the CS contribution in eq. ( |3.36dD by 3/2 one obtains J'^ + J'^ = iVc/(4vr2) Bfi^ and 

- = Nc/{4:iT^) Bfi @, both of which are 1/2 times the results of refs. |] and 
1^], respectively (cf. sec. |0| below). 

In the remainder of the paper we shall again consider the full, unmodified chiral cur- 
rents (^]^) which contain the complete covariant anomaly upon inclusion of the counterterm 

dm. 
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4. Axial and vector currents 

In this section we shall use the results from the previous sections to compute the vector and 
axial currents in the presence of a magnetic field and a quark chemical potential /i as well 
as an axial chemical potential We have seen that a consistent definition of the currents 
is not obvious in the given setup. We shall focus on the definition of the currents presented 
in sec. |2.2| since they reproduce, together with Bardeen's counterterm, the correct anomaly. 
Before going into the details, let us explain the expected physics behind the vector current, 
in other words the chiral magnetic effect. 

4.1 The chiral magnetic effect 

A (noncentral) heavy-ion collision, where the chiral magnetic effect is expected to occur, is 
more complicated than we can capture with our thermodynamic description. The physical 
situation and its simplified description within a thermodynamic approach is as follows 
||2[. In the high-temperature phase gluonic sphaleron configurations with nonzero winding 
number should be produced with relatively high probability, inducing an imbalance in left- 
and right-handed quarks due to the QCD anomaly and thus a nonzero axial number density 
ns . In the simple picture applied in ref. , such chirality changing transitions are assumed 
to have taken place in a nonequilibrium situation, after which in equilibrium a finite is 
no longer changed by the QCD anomaly. The QED anomaly on the other hand does not 
change as long as only a magnetic field is present, so can be considered a conserved 
quantity for which we may introduce //s as the corresponding chemical potential. (We have 
introduced also electric fields above for the sake of checking the axial anomaly, but shall set 
them to zero in the final results.) Nonzero quark masses and/or nonzero chiral condensates 
can be expected to lead to a decay of n^. In the given context, it is thus questionable to 
apply the equilibrium description also to the chirally broken phase, and strictly speaking 
our approach should be extended to a nonequilibrium calculation. 

Let us now briefly recapitulate the physics behind the occurrence of the vector current 
which constitutes the chiral magnetic effect in terms of a (quasi) particle picture |2|. 
Suppose the magnetic field leads to a spin polarization of all fermions, i.e., the spins of all 
quarks are aligned parallel or antiparallel to the magnetic field depending on their charge 
being positive or negative. Massless right-handed fermions, which have positive helicity, 
have momenta parallel to their spin, so they move parallel to the magnetic field if they 
have positive charge, and antiparallel otherwise. For left-handed fermions with negative 
helicity, the situation is exactly reversed. If there are more right-handed than left-handed 
fermions, > 0, there is a resulting net electromagnetic current parallel to the magnetic 
field. (Antifermions have helicity opposite to chirality but also opposite charge, so they give 
a current in the same direction.) For weakly-coupled fermions this picture applies since 
in the lowest Landau level indeed all fermions have their spins aligned in the direction of 
the magnetic field according to their charge. The chiral magnetic effect then results solely 
from the lowest Landau level. The contribution of fermions in higher Landau levels, where 
both parallel and antiparallel spin projections are populated, cancels out. This can be seen 
explicitly upon using the thermodynamic potential of free fermions in a magnetic field, and 
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the resulting current is 



J 



(4.1) 



In our model we cannot see any Landau levels directly. Therefore, let us also repeat another, 
apparently more general, derivation of the chiral magnetic effect. It is based on an energy 
conservation argument originally pointed out by Nielsen and Ninomiya [45| and applied 
in ref. It states that an energy 2^5 is needed to replace a fermion at the left-handed 
Fermi surface /ul with a fermion at the right-handed Fermi surface fiR. This conversion 
changes the axial number density by dN^ = 2, i.e., the energy actually is fi^dN^. Such a 
change in A^5 is possible through the QED anomaly in the presence of an electric and a 
(non-orthogonal) magnetic field. The energy can thus be provided by an electric current. 
Hence, the change in per unit volume and time is given by the electric power per unit 
volume J • E, 



3 E 



dt 



(4.2) 



Now we know from the (covariant) axial anomaly ( p.l9bD that, with V- J5 = and = J^, 



dn^ 
~dt 



Nc_ 
27r2 



BE. 



(4.3) 



Inserting this into eq. ( [4.2| ) and taking B and E parallel yields a current J in the direction 
of B, given by eq. ( |4. 1| ) . Note that this argument only works for a nonzero, although 
arbitrarily small, electric field. 

Besides the vector current we shall also compute the axial current for which the anal- 
ogous topological result is ^ 



J5 



Nc_ 
27r2 



(4.4) 



which is proportional to the ordinary quark chemical potential and thus of potential interest 
in neutron and quark star physics. 

4.2 Currents w^ith consistent anomaly 

We have already computed the currents in the symmetric phase, see eqs. (I^H). The 



analogue for the broken phase, obtained from the definition (2/7) and the gauge fields and 

Et±{jt-ety 



field strengths ( pi^ ) and (|3.7[ ) is 



J. 



KK 



L/R - 3 
3 _ -^/i„)i/r2 



PX2 



//5,t coth Pit =f {jt - et) + 
P coth fdn 



JE/R = T4acM^k/? 



=FA'5,t - (jt - et) coth Pit 



(4.5a) 
(4.5b) 
(4.5c) 

(4.5d) 



Again, to make the origin of the various terms transparent we have written the CS con- 
tributions separately. All terms containing a 1/3 come from the CS action. As for the 
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Figure 3: Vector current J\\ = jTa per imbalance of right- and left-handed fermions = as a 
function of the dimensionless temperature — 2ttT/Mkk for different values of the dimensionless 
magnetic field (3 = aB/M^^ ~ 5/(0.35 GeV^) ~ 5/(2 • lO^^G) (left panel) and as a function of 
/3 for different values of 9 (right panel). The critical temperature for the chiral phase transition is 
Tc = MKK/(27r), i.e., 9c = 1. The currents in this plot are obtained using the consistent anomaly, 
i.e., before adding Bardeen's counterterm to fulfill the covariant anomaly. After this term is added, 
the vector current vanishes exactly. The left plot shows the discontinuity at the first order chiral 
phase transition. This discontinuity vanishes for asymptotically large magnetic fields. The right 
panel shows that the current saturates at a value of J|| = ^J^^, in very good agreement with the 
lattice data for the root mean square value of fluctuations of vector currents and axial densities 
p^ . The three lattice data points are taken from from figs. 4 and 8 of ref. and correspond to 
a temperature T = 1.12 Tc. The shaded area indicates the results read off from fig. 11 of ref. ||2^ 
for the cleaner case of a (T = 0) instanton-like configuration, where the corresponding points lie 
between J\\/J^ — 0.66 — 0.77 for magnetic fields of /3 ~ 3.0 and higher. 



symmetric phase, we can easily check the consistent anomaly ( 2.10| ). Using the expression 
for the supercurrent ( |3.11| ) and = NcB/lGir'^ we find 



and 



Nr Nf. „ /3coth/37r\ , , „, 

^I^^.Tf^I, = ^j^B [e ^ ) . (4.7) 



which confirms eq. ( |2.1[1| ). The axial electric field seems to be modified by a complicated 
function of the dimensionless magnetic field. This originates from the mixing of the electric 
field with the supercurrent, which both enter the boundary value of ^3. We shall see that 
this somewhat strange structure disappears after adding Bardeen's counterterm. 

From the results ( p.36 ) and ( [4.5[ ) we may compute the vector currents in the chirally 



symmetric and broken phases. For the following results we set E = e = 0. We find the 
same result for both phases which is 

Js = (Jym + Jcs)3 = (1 - ^) ^^/'s • (^-^^ 
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This differs by a factor 2/3 from the topological result (4J-). This difference is not surprising 
since we have not implemented the covariant anomaly yet. To this end we must add 
Bardeen's counterterm. Before doing so we point out an interesting result which we obtain 
by considering the ratio of the vector current over the axial density. From eqs. ( |3.36| ) and 
( |4.5| ) we obtain 

2/3 go 



^ ^ 2 
3 



^„ for T > Tc 

6*^ Po , (4.9) 

tanh /Jvr for T < Tc 



which is displayed in fig. |3|. In the left panel we see that the first order chiral phase transition 
manifests itself in the discontinuity of the ratio J^/^J^- Interestingly, the jump vanishes for 
asymptotically large magnetic fields. The curves for the symmetric phase are in qualitative 
agreement with the weak-coupling results in fig. 2 of ref. [0]. The right panel shows an 



intriguing agreement of our result for the ratio Jzj with recent lattice results [22| for the 
root mean square values of electric currents and chiral densities at large magnetic fields. 
While the very good numerical agreement might be a coincidence, the lattice results as 
well as our result clearly show an asymptotic value significantly smaller than 1. If it were 
1, the entire imbalance in right-and left-handed fermions, i.e., all excess right-handed 
fermions, would contribute to the current for asymptotically large magnetic fields. This 
is expected at least at weak coupling. In this case, for sufficiently large magnetic fields, 
all fermions populate the lowest Landau level. Consequently, since the current originates 
solely from the lowest Landau level, as explained above, one expects JzU^ 1. This 
is confirmed in the weak-coupling calculation of ref. |2[, see fig. 6 in this reference. The 
lattice result suggests that at strong coupling there may be important modifications to 
the Landau level picture. We emphasize, however, that fig. ^ is not yet our final physical 
prediction. The model has not yet been appropriately renormalized in order to exhibit the 
covariant anomaly. 

We also remark that the scale of our magnetic fields is very large such that for all phys- 
ical applications, be it in heavy-ion collisions or in magnetars, the limit of weak magnetic 
fields is sufficient. In fact, a dimensionless magnetic field P = 1 corresponds roughly to a 
magnetic field 5 ~ 2 • 10^^ G if one follows refs. ^ and sets Nc = 3, Mrk ^ 949 MeV, 



K ~ 0.00745, which fits the experimental values for the pion decay constant and the rho 
meson mass. Therefore, for all applications we have in mind, /? ^ 1. Moreover, one should 
recall that we have used the YM approximation for the DBI action. This is of course a 
good approximation for small magnetic fields, but our extrapolation to larger magnetic 
fields may be subject to modification when the full DBI action is employed. On the other 
hand, in the limit /? ^ 1, the results for the on-shell action, eqs. ( A.13| ) and ( A.33| ) exhibit 



a strong suppression of the YM action compared to the CS action. This suggests that our 
approximation is reliable also for asymptotically large magnetic fields. 

4.3 Currents vi^ith covariant anomaly and absence of the chiral magnetic effect 



The next step is to include Bardeen's counterterm ( 2.13 ) in order to implement the covarl 



ant anomaly. In the broken phase there is a slight complication because the counterterm 
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4^2^ (3 qo 


4vr^ P go 







T<T, 


Gvr^ p 


^^/is-B coth/Jvr 





Nc ^fi coth /Jtt 



Table 1: Vector and axial densities ^ and vector and axial currents jTji, in the direction 

of the magnetic field B after adding Bardeen's counterterm. All results are given as functions of 
the dimensionless temperature Q = 27rT/AfKK and the dimensionless magnetic field /3 = aS/M^j^. 
The densities in the chirally symmetric phase (T > Tc) depend on temperature; the ratio po/qo 
behaves as po/'Zo 1 for (3/0'^ and po/qo 2(3/9^ for (3/0^ — + cx3. In the chirally broken phase 
(T < Tc), all quantities are independent of temperature; the function p behaves as p — s- 3/{2tt) 
for (3 —^ and p ^ (3 for /3 — > oo. The vector current vanishes exactly in both symmetric and 
broken phases; this indicates the absence of a chiral magnetic effect in the Sakai-Sugimoto model, 
see discussion in the text. For the axial current, the temperature-independent topological result 
is reproduced in the symmetric phase. See fig. ^ for the comparison of the axial currents in the 
symmetric and broken phases. 



should only involve genuine background gauge fields, and not those boundary values of 
the bulk gauge fields that due to the gauge choice = Q represent gradients of the pion 
field. This means that we have to subtract the time-independent part of the supercurrent 
J = —(3p/2p from the boundary values of the A^^^ ■ Then, with eq. ( |2T5| ) and the value 
of c from eq. (|2.18 ) , the contributions of the counterterm to the currents are 



T <Tc: 
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_ 4kM2j, 
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(3 coth Ptt 
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(3 coth fSn 



2p(/3) 
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(4.10a) 
(4.10b) 
(4.10c) 



:^^/3(^,±2/.5,t).(4.10d) 



The first observation is that the 2-component of the current vanishes after adding the coun- 
terterm. As mentioned above, this 2-component was unphysical anyway. The cancellation 
of this component is therefore, besides the covariant anomaly, another sign for the necessity 
of the counterterm. The covariant anomaly is now correctly contained in the renormalized 
currents i7^/^ = c7/^/r + ^•^un- This is clear by construction, and can also be verified 



explicitly: adding eqs. (|4.10D to eqs. ( [4. 5]) , yields 

df,J^ = , d^J^ = 
with the vector and axial currents . 



Nc 



dtJl = ^BE, 



(4.11) 
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Figure 4: Axial current in chirally symmetric (T > Tc) and chirally broken [T < T^) phases in 
the presence of chemical potential fji as functions of the dimensionless magnetic field (3 = aS/M^j^, 
i.e., the magnetic field in units of M^j^/a ~ 2 • 10^^ G. In the symmetric phase the current is 
linear in B, while the current in the broken phase is linear only for small B and asymptotically 
large B as indicated by the dashed lines. Due to the huge scale for the magnetic field, the small-S 
approximation for the axial current is sufficient for astrophysical applications. In this case the 
current in the broken phase is simply 1/3 of the current in the symmetric phase. The analytical 
results are given in table |^ where it is also shown that the vector current vanishes. 

The contributions of the counterterm to the currents in the symmetric phase are 

r > T, : Ajo ^ = ±^^iMk^(3^J/^ _ j^L/R^ ^ ± 26) , (4.12a) 

l^Jl/R = Q, (4.12b) 
^Jl,n = T-^^Px,{E^e), (4.12c) 

AJ,% = ^^/?(3</^ - A'J^) = ^i^/3(^. ± 2^,0 .(4.12d) 

These counterterms have to be added to the currents ( |3.36D . Again, the 2-component of 
the currents is canceled, and the covariant anomaly can again be verified explicitly. 

The results for the currents after setting E = e = Q are given in table |l] for both the 
symmetric and the chirally broken phase. For the axial current we find that the counterterm 
exactly cancels the CS part of the current, 

Jl=jl+^jl = {jyy,)l. (4.13) 

In the chirally symmetric phase, this yields exactly the expected topological result ([4.4|). 
In the broken phase, the current is suppressed (but nonvanishing, in contrast to the results 
obtained with a modified CS action [^]). To lowest order in the magnetic fields as well 
as for asymptotically large magnetic fields this suppression is simply given by a numerical 
factor. For intermediate magnetic fields the difference to the symmetric phase is given by 
a complicated function of B. We plot this result in fig. ^. 
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Table 2: Summary of results for the dilferent (parts of the) chiral currents: the Yang-Mills part 
Jym (exclusively considered in ref. the complete current prior to renormalization J7ym+cs 
(= J in the text), the complete current plus Bardeen's counterterm J7ym+cs + A J' = J, and the 
chiral current obtained by modifying the Chern-Simons action according to ref. ^ , i/ym+cs ■ 
The correct result for the covariant anomaly is underlined. A "1" in the results for the axial 
current J7|^ means agreement with the exact QCD result of ref. Q; a "1" in the results for the 
electromagnetic current J7|| means agreement with the weak-coupling approach of ref. 



The most striking of our results shown in table |l| is that for both phases the renormal- 
ized vector current is zero for all magnetic fields, 

Js = {Jym + Jcs + ^Jh = " ^ " ^) ^^^^ = 0, (4.14) 

i.e., the chiral magnetic effect has completely disappeared after adding Bardeen's coun- 
terterm. The vector current has been computed in the Sakai-Sugimoto model before, and 
both existing (but differing) results are nonvanishing. One of the results |jl8| is 1/2 of 
the topological result ( |4.1| ). This result, however, has been obtained with the modified 
action discussed in sec. ^ which amounts to multiplying the CS contribution by 3/2 (and 
leaving out the counterterm). As we have seen, this modified action leads to a vanishing 
anomaly. Another result has recently been presented in ref. Q as a limit case of a more 
general frequency-dependent calculation, but using only the YM part of the current. This 



gives the topological result (|4.1|) , as can also be seen from eq. ([4.14|) . However, as we have 



shown, this does not produce the complete covariant anomaly, see eqs. ( p. 21 ). 



One of the purposes of our paper is to point out the differences of these results and 
the problems of the various approaches regarding the correct anomaly. A summary of our 
findings is given in table H Although in our approach the correct anomaly is ensured, we 
do not claim to have the final answer since the problem of the ambiguity of the currents, 
see sec. remains. Our approach shows that the CS part of the currents is important 



for two reasons. First, as realized already in earlier works 23, 26 1, it naturally gives a 
nonzero contribution when the currents are defined by varying the full action. Second, and 
maybe more importantly, only by including the CS contribution does one reproduce the 
standard result for the consistent anomaly. And only then one can completely implement 
the covariant anomaly (i.e., a conserved vector current even in the presence of axial gauge 
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fields) by adding an appropriate counterterm as a holographic renormalization. We have 
explained why this counterterm, even in the absence of axial field strengths, but in the 
presence of an axial chemical potential, changes the result for the vector current. We do 
not, however, see a general reason why the counterterm must render the vector current 
zero, i.e., why by requiring the current to be conserved the current itself should disappear 
as it turned out to be the case in our explicit calculation. 

After having understood the difference of our result to previous results in the same 
model, let us discuss its significance in view of the apparent contradiction to the result 
(|4.l| ). As explained above, this result can be derived by using the Landau- level structure 
of fermions in a magnetic field. One might thus view our result as an indication that 
there are no fermionic quasiparticles filling Landau levels in the Sakai-Sugimoto model. 
This may be particularly interesting in view of the recent attempts to see Landau-level-like 



structures in holographic models 47 1. However, as we have pointed out, the derivation 
of the chiral magnetic effect via the energy conservation argument by Nielsen and Ninomiya 
appears to be more general. Obviously, the energy conservation ( [4.2| ) does not hold with 
our results because the left-hand side is zero while the right-hand side yields the expected 
nonzero result from the anomaly, see eq. ( [1.11 ) . More precisely, one can check that eq. 



holds before adding Bardeen's counterterm while the counterterm itself violates eq. ([4.2|). 
However, the form of the counterterm seems to be uniquely determined by the requirements 
of parity and the possibility of accommodating it in holographic renormalization. This 
raises the question whether the apparently general energy argument actually uses properties 
of the system which are different in our strong-coupling approach. Clearly, also in our 
system, chirality is converted by a rate simply given by the anomaly. Possibly the energy 
needed for this conversion cannot be written as in eq. (|4.2| ). A reason might be that this 
energy makes use of the existence of Fermi surfaces for the right- and left-handed particles 
which are absent in our model. It is tempting to speculate that the chiral magnetic effect 
indeed vanishes in the strong-coupling limit and that the weak-coupling results together 
with the recent observations of charge separation in heavy-ion collisions suggest that the 
quark-gluon plasma generated in such a collision is sufficiently weakly coupled to exhibit 
the chiral magnetic effect. A deeper understanding of our result, however, seems required 
before drawing this conclusion. 

We recall that in the context of heavy-ion collisions the magnetic field clearly is time- 
dependent, in contrast to our assumption of a constant magnetic field. Therefore, in order 
to compute the induced current, one has to consider the frequency-dependent chiral con- 
ductivity ^, whereas our result only corresponds to the zero- frequency limit^. In other 
words, even if the conductivity at zero frequency vanishes, a nonzero (time-dependent) 
current can be expected if there is a nonvanishing conductivity at nonzero frequencies. 
However, this does not imply electric charge separation because the separation of charges 
is proportional to the zero-frequency limit of the conductivity H. This is easy to under- 



Judging from the calculation of the chiral magnetic conductivity in ref. (where only the YM part 
was taken into account), one might expect that the full result, taking into account also the (frequency- 
independent) CS part and Bardeen's counterterm, leads to a nonzero conductivity for asymptotically large 
frequency |48| . This seems curious, although we do not see a fundamental reason for this to be unphysical. 
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stand in analogy to a capacitor which cannot be charged with an alternating current, i.e., 
integrating the induced current over time will lead to a vanishing charge separation as long 
as there is no direct current. 



5. Summary and conclusions 

We have studied the strong-coupling behavior of chiral fermions in the presence of a chem- 
ical potential and a background magnetic field in the chirally broken and the chirally sym- 
metric phases. To this end, we have used the Sakai-Sugimoto model which is the model 
that at present comes closest to providing a gravity dual to (large- A'^c) QCD. Our focus was 
the calculation of (topological) axial and vector currents, which are direct consequences of 
the anomalies in the model. 

In particular, we have investigated the chiral magnetic effect, which has been studied 
previously in a weak-coupling approach ||2j and on the lattice [^]. We have pointed out 
that a reliable calculation within the Sakai-Sugimoto model requires a careful discussion 
of the QED anomaly in the model. The standard value of the consistent anomaly arises 
naturally in the model for the most straightforward definition of the current. For this result 
it is crucial to include the contributions from the CS term which are sometimes ignored 
in the literature. The covariant anomaly can then be implemented by adding Bardeen's 
counterterm p^], which is also known to be required in chiral models with a Wess-Zumino- 



Witten term [33|, and we have pointed out that this (finite) counterterm has a form that 
appears consistent with the procedure of holographic renormalization. 

After these general discussions we have solved the equations of motion for the chirally 
broken and the chirally symmetric phases explicitly. In our approximation of the DBI action 
to lowest order in the gauge fields, the solutions are completely analytical. Electric (vector 
and axial) fields parallel to the magnetic field have been considered in order to check the 
anomaly explicitly, but they are not needed and set to zero for our physical (equilibrium) 
results, which only require magnetic fields in the presence of chemical potentials. 

In the presence of a quark chemical potential and a large magnetic field, we have 
calculated the axial current, which may be of interest for astrophysical phenomena such as 



pulsar kicks [10|. In the chirally symmetric phase we have reproduced the known topological 
result 1^, 44 1, while in the chirally broken phase, the current has turned out to be smaller 
but nonvanishing. These results can also be obtained by using only the YM part of the 
current, i.e., in the case of the axial current the CS contribution and the contribution of 
Bardeen's counterterm cancel each other. 

This is different for the vector current. In this case, only the YM contribution yields a 
result in agreement with ref. With the full current, and after adding Bardeen's coun- 
terterm, the vector current becomes zero, indicating the absence of the chiral magnetic 
effect. This is in no obvious contradiction to the result obtained in the weak-coupling 
approach Q, since at weak coupling the chiral magnetic effect is a phenomenon that orig- 
inates solely from fermionic quasiparticles in the lowest Landau level, and it is not clear 
whether this structure persists at strong coupling. 
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In comparison to the result from recent lattice calculations [22| we have pointed out 
an intriguing agreement before adding Bardeen's counterterm, i.e., within the consistent 
anomaly. There the vector current per chirality approaches approximately the value 2/3 
for asymptotically large magnetic field. This is clearly different from the weak-coupling 
approach where this ratio approaches 1. This raises the question whether this asymptotic 
value can distinguish between strong and weak coupling. It also raises the question whether 
the lattice result relates to the consistent, as opposed to the covariant, anomaly. 

Although we have implemented the correct covariant anomaly, also in our approach 
a problem remains. Namely, upon computing the free energy explicitly and then taking 
the derivative with respect to the appropriate source, the currents turn out to be different 
from the straightforward definition via the gauge/gravity correspondence. This somewhat 
disturbing discrepancy can be attributed to boundary terms at spatial infinity. We have 
discussed that a previously suggested fix of this problem by modifying the action [^] seems 
to be not acceptable because it entirely eliminates the axial anomaly from the correspond- 
ingly modified currents. Only the YM part of these currents are anomalous, but those 
suffer from the same thermodynamic inconsistency that this modification was meant to 
fix. Because of these ambiguities, further studies are clearly needed (see also ref. for 
other issues concerning the definition of chiral currents in the Sakai-Sugimoto model). 

Quantitative improvements could be achieved by extending our calculation to the full 
DBI action, though they should be minor for magnetic field strengths of practical interest. 
More critical, but also considerably more difficult, would be the generalization of our ansatz 
to allow for inhomogeneous field configurations and/or inhomogeneous solutions. This 
might be required to resolve the ambiguities in the definition of the chiral currents that we 
have discussed, since those are related to spatial surface contributions in the CS action. 
With our present definition of the chiral currents we have been led to question the very 
existence of the chiral magnetic effect in the strong-coupling regime of the Sakai-Sugimoto 
model (which is gravity dual to large- A^c QCD only in its inaccessible weak coupling limit). 
In this context it would be important to understand whether in the strong-coupling regime 
one has Landau-level-like structures, as conjectured in ref. |I9|. This is interesting also in 
view of recent studies in different gauge/gravity models |46, 47 1. 
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A. Solving the equations of motion 
A.l Chirally broken phase 

In this appendix we solve the equations of motion in the broken phase, eqs. ( |3.4| ). The 



equation of motion for Az ( |3.4c ) is trivially integrated with respect to time t to yield 

kdzAo = -2l3te{z) + kd^Ao , (A.l) 

where we have denoted e{z) = —dtA^ and where we have written the t-independent in- 
tegration constant as kdzAo, to be determined below. Inserting this into eqs. ( ^.4a ) and 
(|3lb| ) yields 



dzikdzAo) = 2pdzA3 + 2ptdze, 
dzikdzAs) = 2f]dzAo - {2(3ft 



2Az) 



k{z) 



(A.2a) 
(A.2b) 



Since the left-hand side of eq. ( A.2a| ) does not depend on t, the right-hand side must be 
independent of t too which implies 



d.A 



-tdze + dzAs , 



(A.3) 



where we have written the t- independent part as dzA^. Consequently, eqs. ( A.2a|) and 
( |A.2b| ) become 



dzikdzAo) = 2pdzAs, 
dzikdzAs) = 2l3dzAo-t 



dzikdzc) 



(A.4a) 
(A.4b) 



Now the square bracket on the right-hand side of eq. ( A. 4b ) must vanish because all other 
terms in the equation do not depend on t. This yields a differential equation for e{z). Since 
e{z) = —dtA^, the boundary conditions for A^ (3^) imply e(iboo) = =p (e — ji), where 
we have decomposed the supercurrent as 



Jt=J + Jit, 



(A.5) 



with J, ji being t- independent. With these boundary conditions the equation for e{z) is 
solved by 



e(z) = E 



cosh(2/? arctan z) 
cosh (3-K 



Ji) 



sinh(2/3 arctan z) 



sinh /Svr 



(A.6) 



To find the solution for Aq and A3 we first conclude from eqs. (A.l) and (A.3), 

cosh(2/3 arctan z) 



Ao{t,z) = Ao{z)+go{t) 

^sinh(2/3 arctan z) 



E- 



A3{t,z) = A3{z)-t 



cosh Pit 

cosh(2/3 arctan z] 



31, 



E 



cosh /Jtt 



(e - Ji 



sinh /Jtt 

sinh(2/3 arctan z) 



sinh [iiT 



(A.7a) 
(A.7b) 



From the z-integration in eq. (|A.lD we have obtained a t-dependent integration constant 
go{t). Such a constant is not permissible in ^3 because of the constraint e{z) = —dtA^. 
Integration constants independent of z and t are included in Aq{z), ^3(2;). We shall fix 
fl'o(i) = Jitcoth/Jvr. This removes the supercurrent from the vector boundary value of 
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AQ(t,z). We cannot at the same time remove the axial field e from this boundary value. 
This becomes clear in hindsight after determining jt from minimization of the free energy. 
Only with the given choice of go{t) this minimization leads to a consistent, i.e., time- 
independent, result for j, j\. It is thus unavoidable for the boundary values of AQ{t,z) to 
become time-dependent, 

Ao(t, 2; = ±00) = /it =F /"5,t , (A.8) 

where we defined 



lj,t = fj, + te coth Pit , /i5^t = fi^ + tE tanh /^tt . 



(A.9) 



We have now reduced the equations of motion ( |A.4| ) to equations for Aq{z), ^3(2;) which 
are simply the gauge fields in the absence of an electric field. These equations can be solved 
in general, 



c 



A f ^\ „ " 2/3arctanz , " „2/3arctanz 

a.q{z} = ao — -r-^e + 777:6 



d_ 



d 



j2/3 arctan z 



(A.lOa) 
(A.lOb) 



with integration constants oq, as, c, and d which are fixed by the boundary conditions 
Aq{^oo) = /i =F /^5; ^3(ico) = TJ- The resulting gauge fields Aq{z), ^3 (-2) are then 
inserted into the gauge fields AQ{t,z) A^{t^z) from eqs. ([A.7|) to obtain the final solution 



sinh(2/3 arctan z) 
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sinh /^TT 

For the free energy we also need the field strengths (times k{z)), 



(A.lla) 



(A.llb) 



kd^Ao = -213 
kd,A3 = -2(5 



cosh(2/3 arctan z) sinh(2/3 arctan z) 



sinh /Jvr 
sinh (2/3 arctan z) 



sinh /37r 



+ (Ji - ei) 



sinh /37r 
cosh(2/3 arctan z) 



sinh /37r 



(A.12a) 
(A.12b) 



As a check, we can perform a parity transformation on the gauge fields. With ^ — > +//, 

^ -Ai5, Jt +Jt, B +B, E -E, e ^ +e and z ^ -z we find Ao{t,z) 
AQ(t,z) and A3{t,z) —A3{t,z), i.e., the fields have the correct behavior under parity 
transformations for each t and z. 
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We can now insert the gauge fields and field strengths into the action (3^) to obtain 



the thermodynamic potential = yS'on-sheii- The YM and CS contributions are 



r^YM 



sinh^ /Jvr 



[{jt - etf - 



' /3coth/?7r 



sinh Ptt 



[{jt - etf - 



t\ , 



(A. 13a) 



(A.lSb) 



where the real time parameter t is treated as an external parameter, unrelated to the 
imaginary time r, whose integration is assumed to just give a factor 1/T. In the YM part 
we have dropped the terms oc B^,E^. This vacuum subtraction can be understood in 
terms of holographic renormalization and follows from the renormalization condition that 
the thermodynamic potential be zero for vanishing chemical potentials; for the explicit 
procedure see ref. [pO|. 



A. 2 Chirally symmetric phase 

Here we solve the equations of motion for the chirally symmetric phase, eqs. ( |3.16| ). For 
notational convenience, let us, in this subsection, denote 



I. 

03 • 



(A.14) 



The time-dependence of the gauge fields is treated analogously to the broken phase. Thus, 
eqs. (3.16c) and ( p.l6b| ) imply 



L/R 
: 



and 



L/R 



-tdzeL/R{z) + dzA. 



■TL/R 



(A.15) 
(A.16) 



where Aq^^, A^^^ are constant in t, and where e^/^j = —dtA^'^. Then, eqs. (|3.16a|) and 
(|3l6b|) read 



L/R 



5.(A;oa,iJ/^) = ±2p'd,A 



L/R 



d.ihd.A^/'') = ±2(3'd,A'-^"'-t 



L/R 



(2/3' 



.2^L/R\ 



ko{z) 



dziksdzCL/R] 



(A.17a) 
(A.17b) 



This is analogous to eq. ( [A.4[ ), the only difference being the two functions ko{z) and k3{z) 
instead of the single function k{z). Again the square bracket in eq. ( A.17b| ) has to vanish. 
This yields a differential equation for eL/ji{z) which is solved as follows. With ei/ji = 
ksdzBi/ji one can rewrite this differential equation as 



(A.18) 
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This equation has the two independent solutions 

p{z) = 2F1 



- 16/3'2 + 1 v^l - 16/3'2 - 1 1 1 



2' l + z2 



1 



^1 - 16/3'2 - 1 v^l - 16/?'2 + 1 3 1 



(A.19a) 
(A.19b) 



Consequently, 6/^/^(2;) = Pl/rP{z) + Ql/r1{z), with constants Pl/r, Ql/Ri a-i^d thus 



eL//?(^) = J2[3^ ^^L/R kodzP + Ql/r kod^q) 



(A.20) 



In the following we need the behavior of the functions p{z)^ Qi^), k^dzP, k^dzq at z = 00 
and z = 0. At z = 00 we have 



p{qo) = -kodzq{oo) = 1 , ^(oo) = kodzp{oo) = . 

At 2; = one finds 

Po = P(0) = - 
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vr 
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and 



(A.21) 

(A.22a) 
(A.22b) 

(A.23) 



kodzp{z ^ 0) = (2/3')% In z , k^z ^ 0) = {2(5' fq^ In z . 

The boundary conditions &lir{.z = cxd) = E^e yield Ql/^j = -(2/3')2(^=Fe)- Inserting this 
constant into eq. ( [AI20D , the result into eqs. ( |A.15| ), ( |A.16D , and integrating the resulting 
equations over z yields the gauge fields 



P, 



L/R 



(2/3 



'J^p{z)-{ETe)q{z) 



+ g^/''{t) + A^'''{z), (A.24a) 



rL/R, 



A^/\t,z) = -t 



P, 



L/R. 



(2/3' 



-2 kodzp{z) - {E=fe) kod^q 



+ 4/^(z) 



(A.24b) 



Here, (t) are time-dependent integration constants from the z integration. We proceed 
by solving eqs. (|A.17 ) for Aq^^, A^^^. Recalling that p{z), q{z) fulfill the differential 
equation ( A.18| ) one easily checks that the functions 

I5'/^(z) = ± 2(3' [Cl/rp{z) + DL/Rqiz)] , 



L/R 



(A.25a) 
(A.25b) 



with integration constants a^^^, a^J^ , C^i^ and Di/fi, are solutions of eqs. ( [Al7|) . One 
now inserts these functions into eqs. ( A.24| ) and determines the integration constants 
as follows. First we recall that all constants except for (t) must not depend on t. 
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This will be used repeatedly in the following. Then we require the boundary condition 
A^^^{t,z = oo) = —t{E =F e) which implies -Dl/h = a^^^- Next, we require regularity of 
A^^^{t, z) at z = 0. With eq. ( [A.23| ) we find that A^^^{t, z ^ 0) diverges logarithmically. 
Requiring the factor in front of the In z term to vanish yields the conditions 
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L/R 



Po 



D 



L/R , 



Po 



(A.26) 



For the temporal component we need to require AQ^^{t,z = 0) = |§ which yields 



L/R _ L/R 



%'"^{^) — ^- With these results the boundary value of ^'^(t, z) becomes 



^/^(t, z = oo) = ^2(i'^[DL/R + t{E T e)] . 
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(A.27) 



This result shows that, as in the broken phase, the boundary values of axial and vector parts 
of Aq necessarily become time-dependent. In other words, in the presence of an electric 
field one cannot fix these boundary values to be time-independent chemical potentials. At 
t = we require Aq (t — 0,z — cxd) — fj, ^ ji^. With these initial values we find 

_Po fJ-T 



D 



L/R 



qo 2(3' 

and the time-dependent chemical potentials become 



A. 



with 



^/^{t,z = oo) 

qo 



/X5 + 20 tE^ . 

Po 



(A.28) 

(A.29) 
(A.30) 



fj. + 2p te — , 
Po 

Collecting all the integration constants, we obtain from eqs. ( [Al2^ ) and (|Al25| ) the final 
solution for the gauge fields, 



Ao^^{t,z) = {fit =F /^5,t) 



p{z) - —q{z) 

qo 



A 



L/R 



{t,z) = -t{E^e)± 



kodzP - — (1 + kodzq) 

qo 



(A.Sla) 
(A.31b) 



Again we can check the behavior of the gauge fields under parity transformations. In 
contrast to the broken phase, we have separate right- and left-handed fields which transform 



as AQ^^{t, z) 



A^^^{t,z), and A^'''{t,z) 



.L/R, 



-Af^^{t,z), as it should be. The field 



strengths become 



kodzA^^^ 



ksd^A^^"^ = ±2/3'(/xt ^ /X5,t 



qo 



p{z) 



Po 



iz) 



(A.32a) 
(A.32b) 



Inserting these results into the action, given by eqs. (|1|) and (l3l^ ), yields the YM and 
CS contributions to the free energy, 



-2Ke^Ml^{f4 + ^2 ,)[/o - (2/3')'/; 



3J , 



{fit + flit) h-{2[3'fh 



Po 

qo 



(A.33a) 

+ 2/3't(/ite + ^5,t^)},(A.33b) 
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where we abbreviated the integrals 



h 







IT i^^odzP - —kodzQ 



°° dz 



p{z) - —q{z) 



(A.34) 
(A.35) 



In the hmit (3^1, the combination Iq — (2/3')^l3 0, so that for very strong magnetic 
fields Qqs ^ ^^YM, as is also the case in the chirally broken phase, see eqs. ( |A.13 ). 
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